arXiv: 1506.0055 lv2 [astro-ph.GA] 1 Nov 2015 


Mon. Not. R. Astron. Soc. 000, 000-000 (2014) 


Printed 3 November 2015 (MN style file v2.2) 


Merging galaxies produce outliers from the Fundamental 
Metallicity Relation 


Asger E. Gr0nnow*^, Kristian Finlator^ and Lise Christensen^ 

^ Dark Cosmology Centre, Niels Bohr Institute, University of Copenhagen, Juliane Maries Vej 30, 2100 Copenhagen, Denmark 


Accepted 2015 May 28. Received 2015 May 28; in original form 2014 September 19 


ABSTRACT 

From a large sample of ss 170, 000 local SDSS galaxies, we find that the Fundamental Metal¬ 
licity Relation (FMR) has an overabundance of outliers, compared to what would be expected 
from a Gaussian distribution of residuals, with significantly lower metallicities than predicted 
from their stellar mass and star formation rate (SFR). This low-metallicity population has 
lower stellar masses, bimodial specific SFRs with enhanced star formation within the aperture 
and smaller half-light radii than the general sample, and is hence a physically distinct popu¬ 
lation. We show that they are consistent with being galaxies that are merging or have recently 
merged with a satellite galaxy. In this scenario, low-metallicity gas flows in from large radii, 
diluting the metallicity of star-forming regions and enhancing the specific SFR until the in¬ 
flowing gas is processed and the metallicity has recovered. We introduce a simple model in 
which mergers with a mass ratio larger than a minimum dilute the central galaxy’s metallicity 
by an amount that is proportional to the stellar mass ratio for a constant time, and show that it 
provides an excellent fit to the distribution of FMR residuals. We find the dilution time-scale 
tober = 1.5681^:^27 Gyr, the average metallicity depression caused by a 1:1 merger to be 
a = 0.2480lQ gQ 2 o minimum mass ratio merger that can be discerned from the 

I pk rxi eyry 

intrinsic Gaussian scatter in the FMR to be ^min = 0.2030lQQQgg (these are statistical errors 
only). From this we derive that the average metallicity depression caused by a merger with 
mass ratio between 1:5 and 1:1 is 0.114 dex. 
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1 INTRODUCTION 


Understanding the way that the metallicities of galaxies depend on 
galactic properties and events is crucial in understanding galaxy 
evolution as metallicity is connected to important galactic pro¬ 
cesses such as inflows, outflows and star formati on. A relation 
betwe en galaxy mass and metallicity was found in lUegueux et al.l 
lll979h and measured with high precision based on SDSS Data Re¬ 
lease 2 as a tight r elation between galaxy stellar mass and gas- 
phase metallicity bv lTremonti et alj J 2 OO 4 I) . This mass-metallicity 
relation (MZR) is thought to arise from galactic winds being more 
efficient at blowing metals out of lower ma ss galaxies owing t o 
their shallower gravitational potenti al wells dP ekel & S ilklll986tl . 
The MZR depends on environment (ICooper et all 200S) and red- 
shift dErb et al.ll200^: iMaiolino et al.ll2008l: iHenrv et al.ll2013h . In 
addition Ellison et alj ( l2008ah found that galaxies with high star 
formation rates (SER) showed systematically lower metallicities 
than galaxies with si milar masses but lower SFRs. This was stud¬ 
ied in more detail in iMannucci et alJ ( 1201 Oh (hereafter MIO) who 
binned stellar masses and metallicities by SFR and found that the 
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MZR anticorrelated with the SFR bin for masses M* < 10^°'® M©. 
In light of this, MIO introduced the Fundamental Metallicity Re¬ 
lation (FMR) between stellar mass, SFR and gas-phase metallic¬ 
ity. They interpreted the SFR dependence as being due to contin¬ 
uous accretion of pristine gas from the intergalactic medium rais- 
ing the SFR while diluti ng the metallicity. Concurrently with MIO 
iLara-Lopez et al.l (1201 Oh also investigated the SFR dependence of 
the MZR and also found a relation between mass, SFR and metal¬ 
licity. These studies of second-parameter dependences of the MZR 
all look at galaxy samples that are selected to differ with the param¬ 
eter of interest whose MZRs can then be compared. A complemen¬ 
tary way of examining further dependences of the MZR or FMR is 
to instead select galaxies with abnormally low or high metallicity 
for their masses (and SFR in the case of the FMR) and then check 
what the metallicity offset correlates with. 

IPeeples et al.l ( l2009l) analysed a sample of 42 metal-poor 
galaxy outliers from the MZR and found that all but two of those 
showed signs of interaction. More systematically, MIO produced a 
histogram of the residuals of the FMR, i.e. the difference between 
the measured metallicity and the metallicity predicted by the FMR 
for each galaxy. These residuals closely followed a Gaussian dis¬ 
tribution with the exception of an extended wing of galaxies with 


































2 A. E. Gr0nnow, K. Finlator and L. Christensen 


lower nuclear metallicities than predicted. There it was speculated 
that the bulk of these low-metallicity galaxies were interacting, but 
this was not examined further. 

That interactions tend to dilute nuclear metallicities 
has been found bo t h in observations o f close galaxy pairs 
(lKewlevrt_ay_|200^ lEllison et alJ l2008bl : iMichel-Dansac et'H] 
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and in simulations of 


mergers 


I Monhjoriet^ _ 20101 ; iRupke et all I 2 OI 0 I : IPerez et ^ 1201 ll : 


iTorrev et alJbOljl with these authors finding systematic offsets 
from the MZR of up to a few tenths of a dex towards lower 
metallicities for merging galaxies. The standard explanation for 
this phenomenon is that the centre of the primary spiral galaxy 
experiences a period of strong inflow of gas from the outskirts of 
the galaxy. This inflow is driven by torques ex erted by stars in bar 
instab ilities created by the tidal interactions jMihos & Herngui^ 
Il996h and it will be metal deficient compared to the nuclear metal- 
licity as spiral galaxies have radial abundance gradients with lower 
metal licities at larger radii dZaritskv et al.|[l9^ : iLuck & Lamberd 
l201lh . The inflow also leads to an increase in star formation which 
causes t he metallicity to even tually recover on a time-scale of a 
few Gyr dMontuori et al.ll20l^ . 

In this paper, we ask whether the population of star-forming 
galaxies whose metallicities are significantly below expectations 
based on the FMR can be readily interpreted within the context of a 
simple model for galaxy mergers that dilute the nuclear metallicity 
and boosts its star formation. In particular: 

(1) Can the low-metallicity tail be successfully modelled as being 
due to mergers? 

(2) Do the galaxies in the low-metallicity tail show complementary 
evidence of being in mergers? 

(3) What does this imply about the impact of mergers on galaxies’ 
gas reservoirs and the time-scale over which mergers have such an 
impact? 

We will show that mergers readily account for the observed 
low-metallicity outliers. This enables us to estimate the time-scale 
and magnitude of merger-induced metallicity dilution in a novel 
way. 


In ^we describe our sample selection. In 0]we fit an FMR 
to our sample and examine the differences between the main sam¬ 
ple and the tail. In 0we review our merger model and in 0we 
describe how we find the best-fitting parameters and their values. 
In ^we discuss uncertainties and assumptions in our model and 
compare our results to observations and hydrodynamical simula¬ 
tions. Finally, we summarize our findings and avenues for future 
work in 0 


2 SAMPLE SELECTION 

We need to establish a large sample of galaxies in or¬ 
der to investigate the tail in the residuals of the FMR. We 
used a subset of the Sloan Digital Sky Survey (SDSS) DR9 
jAhn et alj 1201 2h catalogue by the MPA-JHU groupQ avail¬ 
able at http://www.sdss3.org/dr9/algorithms/galaxy_mpa_jhu.php 
where the techniques used to measure emission lines and de¬ 
rive galaxy parameters are also summarized. Half-light radii were 
adopted from the SDSS-DR9 photometric table “PhotoObjAll”. 
The full catalogue contains 1,843,205 galaxies from which we 


selected 167,086 galaxies according to the criteria which were 
adopted from MIO. 

Only galaxies with redshifts within 0.07 < 2 < 0.30 were 
selected to ensure that the 3 arcsec aperture of the spectroscopic fi¬ 
bre covered a significant part of the galaxies. We also demanded 
that Ha was detected at a signal-to-noise ratio of at least 25. 
This ensures a sufficiently high S/N of the N[II] A6584 line that 
is used in many me tallicity calibrations as well as making BPT 
IIBaldwin et al.lll98ll) diagram classification of the galaxies more 
accurate. We selected only galaxies classified as BPT class 1 (star¬ 
forming) or 2 (low S/N star-forming) filtering out AGNs and com¬ 
posite galaxies. Finally, we filtered out galaxies for which it was 
not possible to measure the [OIII]A5007 line which is necessary to 
determine the metallicity using the 03N2 or R23 calibrations (see 
ED; this amounted to 0.4 per cent of the remaining sample. 

Total stellar masses and SFRs were taken from the MPA- 
JHU catalogue. S t ellar m asses were calculated using the method of 
iKauffmann et al.l ( l2003h . while th e aperture-corrected SFRs were 
estimated through the me thod of Brinchmann et ^ ( l2004h with 
slight modifications from ISalim et al.l l l2007h . The masses were 
multiplied by 1.06 and the SF Rs divided by 1.8 to scale them to a 
Chabrier IMF ( Ichabrieill2003l) . Metallicities were measured within 
the aperture only. 


3 DATA ANALYSIS 


3.1 Fitting an FMR 


We had to use a strong-line method to find the metallicities of 
the galaxies in our sample. These methods are calibrated by fit¬ 
ting the relationship between the ratio of two or more strong emis¬ 
sion lines and metallicities inferred directly from electron tem¬ 
peratures of HII regions. While “direct” metallicities found from 
electron temperature measurements are more accurate, they have 
to be calculated from auroral lines that are very weak and can 
only be detected with sufficient signal-to-n oise in the SDSS spec¬ 
tra by stacking ([Andrews & Martinill2013h . This is a problem as 
we need to be able to measure metallicities for single galaxies as 
we wi ll be studying outlie rs. We chose to use the 03N2 calibra¬ 
tion of [Marino et al.l ( l2013h . which gives 12 -f log(0/H) as a linear 
function of 03N2= log ([OIII]A5007/H,9 x Ha/[NII]A6584), to 


derive metallicities as this recent calibration is based on more ex- 
tensive electron tempe rature data than older calibrations such as 
[Pettini & Pagel (2004). MIO used an average of the N2 and R23 
calibrations of Maiolino et alj ( [20081) but as we are interested only 
in the residuals, the differences in the absolute metallicities de¬ 
rived from different calibrations do not matter and the differences 
in the residual distributions are quite small (see ^6.4t . We fitted 
an FMR of the form introduced in MIO to the data using least 
squares (see ^6.61 for a discussion of alternative parametrizations 
of the mass-SFR-metallicity relation). While MIO only fitted their 
FMR to galaxies with > 10®'we did not include a mass 
cut. In [Mannucci et al.[ ([201 ih . where the FMR was extended down 
to M* « 1O®'®M0, it can be seen that the FMR for low-mass 
galaxies begins to deviate significantly from the extrapolation of 
the FMR of MIO at masses below ~ 10®'®Mq. Only 0.1 per cent 
of the galaxies in our sample has such low masses. Our fit yielded 
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(12 -P log(0/H)FMR) = 8.504 -P 0.169m - 0.034s - 0.110m® 

-P 0.082ms-0.048s®, (1) 
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Figure 1. The black histogram is the residuals of the fitted FMR, the unbroken line is the best-fitting model which has t = 1.57 Gyr, a = 0.248 and 
?min = 0.205 and the dashed line is the best-fitting Gaussian. The grey part of the histogram at large positive offsets is the high-r tail that is not included in 
the model. 


where m = log(M,) — 10 and s = log(SFR) and m and s are in 
units of Mq and M© yr“^, respectively. The standard errors of the 
coefficients are all between 0.1 and 1 percent except for the first 
coefficient which has a much smaller error, and the off-diagonal 
elements of the covariance matrix are much smaller than the diag¬ 
onal elements. The residuals r = (12 + log(0/H)data) — (12 + 
log(0/H)FMR) were then computed for all galaxies. 

We show a histogram of these metallicity residuals in Fig. [T] 
This histogram has 100 equally spaced bins from —0.5 to 0.5 dex. 
58 of the 167,087 galaxies, i.e. about 1 in 3000 galaxies, fall out¬ 
side this range. All but a single one of these outliers (which falls 
just below the range with r = —0.5004 dex) have r > 0.5 dex 
and appear to be mostly dwarfs with log{Mf /Mq) ~ 8 and high 
specific SFRs (SSFRs). We will not examine these high-metallicity 
outliers further and from this point on in the analysis these outliers 
are filtered out. 

This metallicity dispersion is almost Gaussian distributed but 
with a distinct tail towards lower metallicities as was also noted in 
MIO (cf. their fig. 3). We fit a Gaussian function to the metallicity 
dispersion using least squares finding a dispersion of cr = 0.048 
dex and a slight offset from zero of /r = 0.007 dex; this enables us 
to quantitatively define the tail as the bins with r < fj,—2a. This tail 
contains an excess amount of 4.29 per cent of all galaxies relative to 


the number of galaxies that it would contain if it followed the fitted 
Gaussian. Upon close inspection a tail towards high metallicities 
can be seen as well but this tail only contains an excess of about 1 
per cent of all galaxies (where this tail is defined as the bins with 
r > fj, + 2a). 

MIO binned their stellar masses and SFRs in bins of 0.05 dex 
and then fitted their FMR to the median values in each bin. Doing 
this in our case turned out to make very little difference in the fitted 
FMR so we chose not to bin the masses and SFRs. 

3.2 Properties of the galaxies in the tail 

We compare the galaxies in the tail to the general population to es¬ 
tablish whether or not they form a distinct population. The stellar 
masses of the galaxies in the tail are generally lower than the gen¬ 
eral population with a median value that is 0.19 dex below the me¬ 
dian value of all galaxies, both being approximately lognormally 
distributed (see Fig. [2j. This suggests that the tail cannot be an 
artefact of metallicity aperture effects alone (i.e. galaxies that are 
more completely covered by the aperture being found to have lower 
metallicities because of metallicity gradients) because in that case 
the mass distributions should be similar as the stellar masses are 
all estimated for entire galaxies. It also indicates that the tail is not 
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Figure 2. Histograms showing the normalized mass distribution of the over¬ 
all sample (red) and only the galaxies in the tail (green, dashed). The blue, 
low-probability curve is the mass distribution of the galaxies in the tail us¬ 
ing the same normalization as for the overall sample to compare the sample 
sizes. 

predominantly driven by errors in mass because in that case the 
galaxies in the tail should have erroneously high masses as this 
would cause their predicted metallicities to be too high. That the 
galaxies in the tail have lower masses than the general population 
can be interpreted as being due to low-mass galaxies being more 
strongly impacted by mergers because of weaker bulges or higher 
gas fractions. 

The aperture-corrected SSFRs of the galaxies in the tail form 
a bimodial distribution with a peak about 0.5 dex below the mode 
for the general SSFR distribution and another, smaller peak about 
0.6 dex above the mode for the general distribution (see Fig. O. 
Overall, the median value of the SSFRs of the galaxies in the tail 
is slightly below (0.08 dex) the median value for all galaxies. A 
simple interpretation of the high-SSFR peak in the tail is that many 
of those galaxies are experiencing a boost in star formation ow¬ 
ing to interactions. In Fig. |4] we show the distribution of SSFRs 
within the aperture only. These were taken from the MPA-JHU cat¬ 
alogue as well and were derived using similar methods but with¬ 
out the photometric aperture correction. These are more accurate 
than the aperture-corrected SSFRs as the SFRs and masses are 
based only on spectroscopy. In this case, there is no bimodiality 
and the SSFRs of the galaxies in the tail are generally higher than 
the general population with a median value that is 0.5 dex above the 
median for all galaxies. Unlike the SSFRs derived in MIO, where 
SFRs were estimated within an aperture while using global masses, 
these aperture SSFRs are physically meaningful. They show the 
relative amount of star f ormation in the sam e region as where the 
metallicity is measured dSalim et al.ll2014l) . Enhanced SSFR (or, 
equivalently, enhanced SFR compared to control galaxies of sim- 
ilar mass) within the SPSS fibre has been linked to interactions 
llLi et alj|2008l : [Judder et ^l2012l : [Patton et al.ll2013h . The origin 
of the bimodial shape of the aperture-corrected SSFRs in the tail is 
puzzling as none of the other parameters show any bimodiality and 
examining the subset of tail galaxies that are very close to one of 
the peaks reveals only slight differences. As there is also consider¬ 
able overlap between the two SSFR peaks there is no simple way 
of filtering out either peak. A morphological study, that is outside 
the scope of this paper, would be necessary to check whether the 
low-SSFR tail galaxies are typically interacting. 

The metallicities of the galaxies in the tail are as expected sig¬ 
nificantly lower than for the general population the median metal¬ 
licity of the tail being 0.17 dex below the overall median value (see 
Fig. [Hi. One might worry that a significant part of the tail would 
then have metallicities below 12 -|- log(0/H) = 8.2, which is the 


lower limit of validity for the 03N2 calibration we use. However, 
only 4.4 per cent of the galaxies in the tail have metallicities below 
this. For the overall sample 99.2 per cent are within the range of 
validity of 8.2 <12-1- log(0/H) < 8.8. 

As can be seen in Fig. the galaxies in the tail have signifi¬ 
cantly smaller half-light radii than the general population. The fig¬ 
ure shows the U-band half-light radii but the G and R bands show 
the same trend. The difference in the median values is about 30 per 
cent which is small enough that it could be attributed to the differ¬ 
ences in the mass and redshift distributions of the tail and general 
population samples. However, the sudden rise at 0.5 arcsec suggests 
that the half-light radii of many of these galaxies are overestimated 
because they are smaller than the seeing. Additionally, as can be 
seen from the figure, the modes of the two half-light radius distri¬ 
butions differ by a factor above 2. Thus a significant part of the 
galaxies in the tail are significantly more co mpact than is typica l 
for the full sample. This is in agreement with lEllison et alj (l2008bl) 
who found that galaxies in their pair sample with half-light radii 
below 3 kpc tended to have lower metallicities than galaxies with 
larger half-light radii. 

In summary, galaxies whose metallicity is lower than pre¬ 
dicted for their stellar mass and SFR have systematically low stellar 
mass, high SSFR within the aperture but slightly lower SSFR when 
aperture-corrected, small half-light radius and low metallicity. The 
tail towards lower metallicities is much larger than the tail towards 
higher metallicities. These considerations support the view that the 
low-metallicity tail consists of a physically distinct population. 


4 THE MERGER MODEL 

We assume that the residuals of the FMR in the absence of merg¬ 
ers would be a normal distribution and fit a Gaussian function to 
the residuals from the FMR found in 43. II using least squares. A 
Gaussian G(r) is then fitted again but this time only in the inter¬ 
val [/i — 2 (t, -b 2a], where the mean /r and the standard devi¬ 
ation a is estimated from the first fit, in order to avoid the low-r 
tail where the shape deviates significantly from a Gaussian and the 
high-r tail which is not a part of the model. We then introduce a 
simple model to take the effect of interaction-triggered metallicity 
dilution into account. This model has three free parameters, ^min, 
a and r, and two fixed parameters, namely the scatter and mean 
of the non-merging population, G(r), that has already been fitted, 
^min represents the minimum mass ratio merger that can be dis¬ 
cerned from non-merging population while a and r represent the 
magnitude and time-scale of the dilution, respectively. 

Mergers above a certain mass ratio threshold, ^min, shift galax¬ 
ies towards lower r by an amount a^*, where is the stellar mass 
ratio of the merger with respect to the most massive galaxy. So in 
the case of a 1:1 merger the metallicity changes by —a dex assum¬ 
ing that the change in r caused by changes in stellar mass and/or 
SFR is small compared to the change in r caused by the metallic¬ 
ity change (we show that this assumption is quite accurate in 46.11 1. 
The metallicity change might of course vary with mass ratio in a 
more complicated way but as we have no knowledge of this rela¬ 
tionship outside of that larger mass ratios should cause more dilu¬ 
tion we assume this simple first-order relation. This offset remains 
constant for a time r (in Gyr) before the galaxy’s metallicity returns 
to normal, i.e. goes back to follow the Gaussian scatter around the 
FMR. ^min can be thought of as representing the minimum shift in 
r that can be distinguished from the intrinsic Gaussian scatter of 
the FMR, so we would expect a^niin ~ a. Note that our definition 
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Figure 3. Same as Fig.[2]but for the aperture-corrected SSFR (SFR/M*). 



log{0/H) + 12 


Figure 5. Same as Fig.[2]but for metallicity. 




Figure 4. Same as Fig.[2]but for the SSFR within the fibre. Figure 6. Same as Fig.[2|but for U-band half-light radii. 


of the stellar mass ratio, £+ = entails that ^ 1, so 

*,primary 

secondary merger members (i.e. the less massive galaxy in the pair) 
are excluded from the model. This is intentional as we do not ex¬ 
pect secondary merger members to experience notable metallicity 
dilution (see ^6.2b . 

The probability density function (PDF) of FMR residuals is 
given by the following expression: 

P{r) = P(r|unmerged)P(unmerged)-|- 

P(r|merged)P(merged). (2) 

Obviously, P(unmerged) = 1 — P(merged). P(merged) is calcu¬ 
lated hy integrating the merger rate per galaxy per mass ratio per 
lookback time for a galaxy of stellar mass M,, ^ P(meiged^g.,M.) ^ 
over lookback times 0 to T 3 > t, mass ratios from to 1 (as 
is defined as the stellar mass ratio of the merger with respect to the 
most massive galaxy) and the galaxy stellar mass PDF . That 
is. 


P( merged) = 



d^P(merged|^*, M,) 
d^* dt 


^ df d^* dM*. 


(3) 


P(r|unmerged) is just the normal distribution previously fitted to 
the residuals, G(r), normalized to have unit area. 

We calculate the second term in P(r) hy integrating the proh- 
abilities given a specific mass ratio over : 


P(r |merged)P(merged) = 


J P(r|merged, ^*)P(merged, ^*) d^*. (4) 


P(r|merged,^*) is just the Gaussian P(r|unmerged) where the 


mean is shifted from /r to /r — a^*. As per equation 0, equation 
@ becomes 


P(r |merged)P(merged) = 


Jo aivu J. Jo 


(5) 

Q^* clc J. 


We obtain the merger rate per galaxy per stellar mass ratio 
per lookback time from the merger rate per halo per halo mass 
ratio per redshift (fo r a halo mass Mh) z) from 

iFakhouri et al.1 fioi^ through the following calculation: 


dV 

d^* df 




d^h dz 
d^P 
d^;, dz 
M. 




d^h dM*,s dz 


iMh,^h,t{z)) 


dM*,s d^« df 
dMhiM,,. 


dM* 


M. 


Mh{M, 


(1 -f z)Ho\/^m -f Da, 


(6) 


where M*,s is the stellar mass of the smaller galaxy in the merger 
and Ho is in units of Gyr“ . We use a ACDM cosmology with 
Ho ~ 70kms“^ Mpc~^, Da = 0.7 and Dm = 0.3 and the stel¬ 
lar mass - halo mass relation of ISehroozi et alj ( |2010|) . Normally 
some “merger delay” time-scale would he imposed because a halo 
will spend some time as a subhalo before the merger is completed. 
We do not include this effect as the metallicity depression may well 
begin and end at different times than the merger itself as estimated 
from e.g. dynamical friction and operate on a time-scale that scales 
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differently with parameters such as mass ratio from the merger de¬ 
lay time-scale (see 36.21 1. 

We evaluate equation da numerically creating a histogram 
of galaxies in bins of width 0.005 dex in r. We do the time in¬ 
tegration by simply multiplying by r under the assumption that 
the merger rate is approximately constant over T ^ t (i.e. as¬ 
suming T tn where tn is the Hubble time) with the look- 
back time set to t{z = 0.1) = 1.30 Gyr, the typical redshift of 
the galaxies in our sample. The stellar mass integration is done by 
looping over 10 stellar masses equally logarithmically spaced from 
log{Mt/M q) = 9.125 to log(M«/M 0 ) = 11.375 with 2^ be¬ 
ing based on the mass distribution of our sample. Finally, the prob¬ 
ability is converted to galaxy counts in each bin in order to compare 
the model to the data. 


5 MODEL FITS 

We fit our five-parameter model in two steps: First, we fit a Gaus¬ 
sian to the FMR residual distribution in the way described in 0 We 
found that the best-fitting Gaussian had a mean of 0.006 dex and a 
scatter of 0.046 dex. Then we fitted the three parameters describ¬ 
ing the low-r tail: ^min which represents the minimum mass ratio 
merger that can be discerned, a which represents the metallicity 
change caused by a merger and r which represents the time-scale 
of the dilution (see 0 for a detailed description of these parame¬ 
ters). We ran 40^ models with 40 equally spaced values of each of 
these. 

We evaluate the relative goodness-of-fit of each model by 
defining a likelihood that compares the predicted and observed 
number of galaxies in each bin of r. In particular, we calculate 
each model’s negative log-likelihood assuming that the errors are 
Poisson distributed. This is a more appropriate goodness-of-fit in¬ 
dicator than in this case because several of the bins in the low-r 
tail have fewer than five galaxies and therefore have errors that are 
distributed significantly differently from a Gaussian and rebinning 
to ensure that all bins had at least five galaxies produced a bin¬ 
ning that was too coarse. The negative log-likelihood for Poisson 
distributed errors is given by 

— lnL= In (di!)-I-rrii — di In mi. (7) 

iG.r bin no. 

Here di is the number of galaxies in bin i and nii is the number of 
galaxies in bin i predicted by the model. The set of parameters r, 
a and ^min that produces the model that fits the data the best are the 
ones that minimize — In L as this is equivalent to maximizing the 
likelihood L (we call this maximum likelihood Lq). As In (di!) be¬ 
comes enormous and is constant across all the models anyway we 
set this term to zero and find the model that minimizes this shifted 
log-likelihood, which will be the same as the one that minimizes 
— In L, instead. We subtract this minimum shifted log-likelihood 
from all the shifted log-likelihoods to get the relative log-likelihood 
A In L = — In L — ( — In Lq) = In 7/0 — In 7/. This is proportional 
to a likelihood ratio test, but we will stick to calling it relative log- 
likelihood and labelling this A In L as likelihood ratio tests are usu¬ 
ally associated with the comparison of nested models. 


5.1 Parameters 

Slices from the AlnL space are plotted in Fig. |7] for r = 
[1.40, 1.712], a = [0.2360, 0.2594] and = [0.1400, 0.2765] 
with 40 equally spaced values in each range, i.e. Ar = 0.008 


Gyr, Aa = 0.0006 and A^min = 0.0035. These parameter ranges 
were found from a coarser initial fitting of models using the ranges 
T = [0.5,4.98], a = [0,0.4985] and Cmin = [0.1,0.88], also 
with 40 equally spaced values in each range. The best-fitting model 
has r = 1.568to.o27 Gyh a = 0.24801^020 and = 
0.2030lo 0095 and is plotted in Fig. □ The given errors are sta¬ 
tistical errors only, the errors owing to various assumptions and un¬ 
certainties in the model are much greater (see As can be seen, 
there are slight degeneracies between r and ^min and between r and 
q; if T is increased a higher ^min and/or a lower a is preferred. This 
behaviour is quite intuitive: a higher number of galaxies showing 
diluted metallicities at a given time because the time-scale of metal¬ 
licity depression is longer can be partially counteracted by exclud¬ 
ing more merger mass ratios or making the effect of mergers on 
metallicity less severe. 

The standard errors are found by marginalising the likelihood 
over the other two parameters for each parameter and finding the 
value of the parameter left and right of the peak where L = 

(Lo being the maximum likelihood) using linear interpolation be¬ 
tween data points. This is done by numerically integrating the rela¬ 
tive likelihood as found from A In L over the other two parameters. 
Using A In L rather than — In L only changes the normalization of 
the marginalized likelihood distribution. For example, in the case 
of a the marginalized likelihood distribution is given by 


L{a) = / / exp(AlnL(a, r, ^mm))drd^„ 

J T J 


// 

J T J £, 


exp(lnL(a,T, ^mm) - InLo) drd^n. 


/ I/(a, T,^,nin)drd^,v 

^0 Jt 


(8) 


can be found by requiring that the integral over all three param¬ 
eters be unity. 

The normalized parameter likelihood distributions are shown 
in Fig.[ 8 | As can be seen, the distributions are approximately Gaus¬ 
sian and quite well resolved, meaning that the simple method used 
to derive the standard errors is appropriate, an assertion that is also 
supported by the fact that integrating from the left or right until 15.9 
per cent of the area is enclosed yields nearly the same estimate of 
the errors. 

To get a feel for how each of the parameters affects the model, 
we change one of the parameters by ±50 per cent while keeping 
the other two parameters at their best-fitting value and the normal¬ 
ization fixed. The resulting models are plotted in Fig. [^ As can 
be seen, increasing r amplifies the tail while decreasing it brings 
the distribution closer to the best-fitting Gaussian. This is because 
increasing r means that galaxies have diluted metallicities for a 
longer time and therefore the relative number of galaxies showing 
depressed metallicity at any one point in time increases. Varying 
a changes the overall shape of the tail with larger values caus¬ 
ing a decrease of the moderate low-r part and an increase of the 
more extreme offsets and smaller values bringing the distribution 
closer to a Gaussian that is wider than the best-fitting one. This 
is because increasing a increases the magnitude of metallicity di¬ 
lution moving galaxies that are already experiencing dilution fur¬ 
ther towards lower metallicities. Increasing/decreasing ^min dimin¬ 
ishes/amplifies the moderate part of the tail while having no influ¬ 
ence at the more extreme offsets. This is because increasing ^min 
removes the galaxies at the lowest mass ratios previously included 
which therefore have the smallest metallicity dilutions. 
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Figure 7. Contour plots of the A In L parameter space with one parameter set to its best-fitting value. The crosses mark the minimum while the black 
overplotted contours mark l-cr, 2-cr and 3-a deviations from the minimum (i.e. A In L = 0.5, A In L = 2 and A In L = 4.5 respectively). 



ct T f ■ 
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Figure 8. Marginalized parameter space for a, t and from left to right. The unbroken line marks the parameter value that yields the maximum likelihood 
Lq (with no interpolation between the parameter values used in the simulation). Dashed lines mark the formal standard errors as estimated by the linearly 
interpolated parameter values left and right of Lq that lead to L = e^/^Lo- 

6 DISCUSSION 

6.1 Magnitude and time-scale of the metallicity dilution 

The magnitude of the nuclear metallicity dilution caused by a 
merger is connected to the parameter a in our model. Because the 
FMR depends on mass and SFR, the change in r in our model is 
affected by changes in mass and SFR occurring during a merger, 
in addition to changes in nuclear metallicity. In practice, however, 
the dependence on mass and SFR of our fitted FMR is sufficiently 
weak compared to the differences in the tail relative to the gen¬ 


eral population that we can safely ignore this. The difference in 
median SFR between the general population and the tail by itself 
only leads to Ar = 0.004 and the difference in mass only leads 
to Ar = —0.026. This is large compared to the statistical errors 
on a quoted in 35.11 but insignificant compared to the systematic 
uncertainties discussed later in this chapter. Thus, we can equate 
the change in r with the change in nuclear metallicity without in¬ 
troducing any significant uncertainty. 

In addition, we lo ok at the sample of 42 metal-poor galaxies 
of IPeeples et al.l (l2009l). These galaxies were selected to be low- 
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12+log(0/H) offset 12+log(0/H) offset 12+log(0/H) offset 

Figure 9. The best-fitting model (unbroken, red line) has one of its parameters decreased by 50 per cent (dashed, blue line) and increased by 50 per cent 
(dot-dashed, green line) while the two other parameters and the normalization are held fixed. From left to right r, a and are varied. 


metallicity outliers from both the luminosity-metallicity and mass- 
metallicity relations and 40 of them were found to have disturbed 
morphologies indicative of interactions. 38 of these galaxies are in 
our sample and of those 30 of them are in the low-metallicity tail 
(with most of the rest being very close to the 2cr cutoff). This shows 
that when taking the SFR dependence into account these interacting 
low-metallicity galaxies remain outliers. 

The best-fitting value of a = 0.248 implies that the nuclear 
metallicity will on average decrease by abo ut 0.25 dex during a 1:1 
me rger in agreement w ith the simulations of iMontuori et aP <20101) 
and lRupke et al.l ( 1201 Oh who found metallicity depressions to be in 
the range 0.2 - 0.3 dex and 0.1- 0.3 dex, respectively. Both of these 
were smooth ed-particle hydrodynam ics simulations of equal mass 
mergers with iMontuori et alJ jioiol) being the more sophisticated 
of the two by including star formation and chemical enrichment 
from supernovae. 

We calculate the average metallicity decrease for 
all mass ratio mergers included in a given model as 
/°(,® 5 rP(r|merged)P(merged)c„ormdr where c„orm = 

P(r|merged)P(merged) dr) (see 0 for the defi¬ 
nitions of P(r I merged) and P (merged)) and find that this is 0.114 
dex for our best-fitting model. This is a somewhat more modest 
decrease than the 0.17 dex metallicity difference we find in the 
median metallicities between our general sample (which is very 
close to the difference between the averages) and the galaxies in 
the tail. However, this value is quite sensitive to the definition of 
where the tail begins, which we chose to be 2a below the mean 
r to be able to clearly distinguish the population in the tail from 
the general population. If, for example, the tail is defined to begin 
at 1.5a below the mean r instead, the galaxies in the tail have an 
average metallicity deficiency of 0.14 dex. 

Comparing our average metallicity diluti on of 0.114 dex to ob- 
servations of close galaxy pairs in the SDSS. lKewlev et ^ l l2006l) 
found that pairs with separations < 20 kpc followed a luminosity- 
metallicity (LZ) relation that was systematically shift ed relative to 
the LZ relation for field galaxies by —0.2 dex but lEllison et al.l 
( l2008bt) argued that roughly half of this shift is due to luminosity so 
correcting for this puts their o bserved metallicity di l utions close to 
our value. However, studies bv lEllison et alJ | [2008bh . l^udder et al.l 


( l2012h and iMichel-Dansac et ^ j2008l) find a smaller change of 
—0.03 to —0.05 dex. Much of this discrepancy is due to the 
use of different lowe r mass ratio limits. lEllison et al.l ll2008bl) and 
IScudder et al.l (1201 2l) both included mergers down to a mass ratio 
of 1:10 while our best-fitting model has ^min = 0.205, i.e. only in¬ 
cludes mergers down to a mass ratio of about 1:5. Keeping the other 
parameters at their best-fitting value but settin g ^min = 0.1 yields an 
average metallicity depression of 0.087 dex. iMichel-Dansac et aP 
( I 2 OO 8 I) does not employ a sharp mass ratio cutoff but they do plot 
the MZR for all galaxies versus interacting galaxies for < 0.2 
and > 0.2 with the larger interactions showing much stronger 
dilution (see their fig. 2). 


The parameter r should be the average time from the mo¬ 
ment when the metallicity becomes diluted due to a merger until it 
completely recovers. Unfortunately most numerical simulations of 
merger-induced metallicity dilution are not run sufficiently long to 
estimate this time-scale and as such only yield lower limits. How¬ 
ever we find our best-fitting value of 1.57 Gyr to be in good agree- 
m ent with the typical me tallicity depression time of 2 Gyr found 


Montuori et alj i2()10t) and above the lower limit of 1 Gyr found 
Torrev et al. I 2 OI 2 I). 


In our model, the metallicity depression is constant in time and 
lasts the same for all mergers. However, simulations show that for 
any particular merger both the magnitude and length of metallicity 
depression depends in complicated ways on the orientations of the 
galaxies and on whether they are on retrograde or prograde orbits. 
The metallicity depression time-s cale probably also depends on 
mass ratio (see g6.5t . Furthermore iTorrev et al.l ( 1201 2l) found that 
the initial gas fractions of the galaxies are important with higher 
gas fractions leading to less dilution. A characteristic feature found 
to some degree in all merger simulations is a “double-dip” shape of 
the metallicity as function of time associated with the first and sec¬ 
ond pericentric passages. We chose to ignore these complications 
as their effects are not well known and including any of them would 
add significant complexity to our model ruining its appealing sim¬ 
plicity. Furthermore, because we are using a rather large sample, 
our model should still yield representative average values of the 
dilution magnitude and time-scale. 
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6.2 Merger rate 

Galaxy-galaxy merger rates have significant uncertainties pertain¬ 
ing to the halo-halo merger rate, the stellar mass-halo mass rela¬ 
tion used to convert halo masses to stellar masses and the method 
used to follow sub-haloes or the assumed merger delay, either of 
which are used to conv ert the halo-halo merg er rate to a galaxy- 
galaxy merger rate (see iHonkins et alj|2010bl for a detailed anal¬ 
ysis of each source of uncertainty in merger rates). For a halo- 
halo merger rate der i ved fr om a simulation such as the one we use 
from lFakhouri et all (1201^ (which is based on the Millennium and 
Millennium-II simulations) there are uncertainties from the defini¬ 
tion of mass ratios, the construction of merger trees and the time 
resolution in the simulation. These yield a combined uncertainty of 
factor ~ 2 . 

A source of uncertainty that we introduce in converting this to 
a galaxy-galaxy merger rate is that we do not include any merger 
delay, i.e. we assume that the time that passes from one of the 
haloes to become a subhalo to the merger is completed is the same 
for all galaxies. This merger delay is often calculated based on a 
model of inspiral due to dynamical friction in which case it depends 
on mass ratio, virial radius of the primary galaxy and the energy 
and an gu lar momentum of t he orbit (see e.g. lBovlan-Kolchin et al.l 
( l2008h or ijiang et al.l j2008l) ). Other methods based on characteris¬ 
tic time-scales fo r gravitational or angul ar momentum capture can 
also be used (see llToDkins et al.l 1 1201 Pall ). In general ignoring the 
merger delay time yields fewer major mergers. The differences in 
merger rates arising from using differe nt merger delays or sub halo¬ 
following methods w as examined in iHopkins et alJ (1201 Obh and 
iHonkins et al.l ( l2010al) where it was found that the merger rate de¬ 
rived using no delay or method to follow subhaloes lay within the 
range of merger rates derived using these methods. More impor¬ 
tantly we are interested in galaxies with diluted metallicities due to 
mergers which may well be a process that scales differently with 
various parameters and begins and ends at different times than the 
merger itself as measured from some dynamical friction or group 
capture time-scale. 

At low stellar masses of \og{M,/M q) < 10, the gas fraction 
is significantly greater than at high stellar masses so our using the 
stellar mass ratio underestimates the actual metallicity impact, 
which should depend on the total tightly bound mass, o f low-M* 
mergers by a factor of at least ~ 3 dHopkins et akll^lObh . 

Qualitatively, the absence of red (i.e. passive) galaxies from 
our sample suppresses the merger fraction, as red galaxies have 
a higher merger rate at the me dian redshift of our sample than 
blue (i.e. star forming) galaxies dkin et al.ll2008h . Thus, the actual 
merger rate should be less than what we obtain from equation ID 
and consequently r ought to be larger. However, the fraction of red 
galaxies in the relevant mass range should not be very large, partic¬ 
ularly given that the galaxies in the tail tend to have lower masses 
than in general, so we expect this effect to be subdominant to the 
others evaluated in this chapter. 

A consequence of defining mass ratio as secondary ^ ^ 

is that we are ignoring the secondary members of mergers in our 
model. The degree of uncertainty caused by this depends on two 
factors: the number of such galaxies in our sample and how the 
metallicity of smaller companions in a merger changes. For a given 
stellar mass we can estimate the fraction of mergers as the primary 
galaxy from the merger rate and the galaxy stellar mass function. 
Based on this calculation and the stellar mass distribution of the 
galaxies in our sample we roughly estimate that if the secondary 
galaxy was affected similarly to the primary galaxy r would be 


overestimated by a factor of ~ 3 - 4 in our model. This calculation 
is detailed in Appendix IaI 

The effect on the metallicity of a galaxy caused by it merging 
with a more massive galaxy has not been examined in any detail but 
we would expect enrichment rather than dilution in this case as the 
smaller galaxy accretes more enriched gas from the larger galaxy 
which also triggers star formation. Indeed IScudder et al.l l l 20 I 2 h 
found that both members of m erging pairs experienc e d simi lar lev¬ 
els of SFR enhancement and iMichel-Dansac et'H] ( l2008h found 
that while massive galaxies in mergers showed diluted metallicities 
smaller merging galaxies with log(M*)/MQ 10 showed enrich¬ 
ment instead. Therefore we expect the vast majority of galaxies that 
are merging or have recently merged as a secondary member to not 
be in the low-r tail and the effect on our results of ignoring these 
should be small. 

We can roughly estimate the merger fraction of our sample as¬ 
suming that all galaxies in excess of the best-fitting Gaussian for 
the FMR residuals below the mean are merging. We find a merger 
fraction of 6.0 per cent of the sample in this way. As we are ignor¬ 
ing flybys (see 36.3t and metallicities recover on a time-scale that 
is longer than the time-scale of the merger being visible (as either 
clear morphological disturbances or galaxy pairs) this is an upper 
limit. We cannot compare this directly to observational measure¬ 
ments of the merger fraction as these are measured for high masses 
(typically > 10 ^° M©) or luminosities only and for different 
types of samples (as we Altered out AGNs and galaxies without ac¬ 
tive star formation) at typically significantly higher redshifts. How¬ 
ever, our value lies within the range of results in the literature found 
from morphologica l merger indicato rs (generally 1-10 per cent at 
lower redshifts; see lLotz et al.l l l201lh ). 


6.3 Flybys 

A potentially important complication that we are ignoring is fly¬ 
bys, i.e. when a galaxy passes another galaxy at a small distance 
but does not merge. We chose not to include flybys as the cur¬ 
rent quantitative knowledge of the effect and rate of flybys is poor 
and because of the considerable complexity it would add to our 
model. Little work has been done to examine the effects of flybys 
on metallicity. However iMontuori et alj 1 I 2 OIOI 1 found that in their 
simulations of equal-mass interactions close flybys caused almost 
as much nuclear metallicity dilution as mergers and that this dilu¬ 
tion lasted for almost as long as for mergers. 

The flyby rate has a lso r eceived very little study. 
ISinha & Hollev-Bockelmai^ ( 1201 2h found from cosmologi¬ 
cal N-body simulations that the rate of “grazing” flybys (i.e. when 
two primary haloes approach each other, overlap for at least half a 
crossing time, then continue on different trajectories as two distinct 
primary haloes again) was comparable to the merger rate for halo 
masses log(Mj,/M 0 ) >11 (corresponding to log(M«/M 0 ) > 9) 
at z ^ 2 . 

In the simple scenario where the flyby rate is equal to the 
merger rate and flybys dilute the nuclear metallicity in a way sim¬ 
ilar to mergers and on the same time-scale, the full effect of ig¬ 
noring flybys in our model is simply that r will be overestimated 
by a factor of 2. While this simple situatio n might not be t hat far 
fro m the truth considering the findings of IMontuori et al.l ( l 2010 h 
and lSinha & Holley-Bockelmannj ( 1201 2h we expect the magnitude 
of dilution to depend on the pericentric distance of the flyby with 
greater distances causing less dilution. 























































10 


A. E. Gr0nnow, K. Finlator and L. Christensen 


Calibration 

Aperture-corrected SFR 

r (Gyr) 

ol (dex) 

'Cmin 

03N2 

Yes 

1 |ry-|-0.029 
J--0/ _0.027 

n o^Q+o ooi’^ 

U.^4«_q QQ20 

u.zuo_0 0115 


No 

-, 7Q+0.046 
-*-• <y_0.042 

u.zuyo_Q 0019 

r\ QQ-|-U.0i49 

y-o^_0.0135 

R23 

Yes 

1 qc-I-0.038 

-L-yo_o.o4o 

0 350+0.0033 
U.OOU_o 0025 

n 909+0.0144 
u.zoz_0 0098 

No 

q'7“I"0.035 

J--O'_0.043 

0 345+0-0023 
u.o^o_o 0034 

U.Z9U_o oil 

N2 

Yes 

-| 11 o-|-0.023 

l.ilO_o.o21 

U.z^oo_o.ooi9 

u.z^o_o.oi36 

No 

-L-^^_0.034 

u.z±u_o 0025 

n 009+0.0189 
u.ooz_o 0159 

N2-tR23 

Yes 

J--oy_o.oi7 

n 9QQ9+0-0017 

U.z»y4_o 0020 

n i 9’7+0.0060 
^•l^^_0.0077 


No 

-L-^^_0.034 

u.z±u_o 0025 

n oo9”t"U.ul?5y 

u.ooz_0 0159 


Table 1. Best-fitting parameters for models based on data using different metallicity calibrations and with and without aperture correction of SFRs. 



12-l-log(0/H) offset 


Figure 10. The FMR residuals of different metallicity calibrations. The unbroken black line is the residuals of the 03N2 calihration of [Marino et alj |2013h 
that we have been us ing while the dashed red, dot-dashed green and dotted blue lines are residuals of calibrations based on N2 iDenicol6Trani^Q2l) . R23 
iMaiolino et al.l2do3l and the average of those two, respectively. 


6.4 Metallicity calibrations 

As different metall icity calibrations y i eld qu ite different metallici- 
ties and MZRs (see lKewlev & EllisonI l l2008i) l we expect the choice 
of calibration to have an impact on the FMR. To evaluate the im¬ 
pact of different calibrations we compare the shape of the distribu¬ 
tion of residuals of the fitted FMR and nature of the low-r ta il of 
the data based on the N2 calibration of lPenicolo et alj l l2002ll . the 


R23 calibration of iMaiolino et al.l l l2008l) and the average of these 
two which was use d in MIO with the p revious data based on the 
03N2 calibration of iMarino et alJfcOl^ . We show the FMR resid¬ 
ual distributions for these calibrations in Fig.llOland the best-fitting 
parameters of the corresponding models in Table[T] 

As with 03N2 these calibrations have limited ranges of valid- 
ity, though both are broad. The range given for N2 i mPenicolo et all 
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ll2002h is 7.2 < 12 + log(0/H) < 9.1 and the entire sample 
lies within this. Recently, however. iMorales-Luis et al] ( l2014h have 
shown that N2 does not correlate with metallicity for galaxies with 
12 + log(0/H) < 7.6. This is not an issue for our sample as only 
three galaxies have such low metallicities. The range for R23 is 
about 7.0 < 12 + log(0/H) < 9.3 and 99 per cent of the sample 
lies within this range. 

iDenicolo et all ( l2Q02h find a relation for 12 + log(0/H) as a 
linear function of N2 = [Nil] A6584/Ha. Using this relation to 
find the metallicities and fitting an FMR to these we find that the 
metallicity dispersion is not too dissimilar to the dispersion based 
on 03N2. Fitting a Gaussian distribution with least squares as be¬ 
fore we find a slightly lower scatter and mean of a = 0.042 dex 
and a fj, — 0.004 dex, respectively. The tail is smaller though, con¬ 
taining an excess of about 3.4 per cent of all galaxies. There is a 
systematic offset towards lower metallicities of 0.05 dex for the 
entire population and the tail alike. The SSFRs of the galaxies in 
the tail have the same bimodial distribution as for 03N2 but with 
a slightly lower median SSFR (0.03 dex lower than for 03N2, i.e. 
0.11 dex lower than the median SSFR of the general population). 
As can be seen in Table [T] the best-fitting value of a is almost the 
same as for 03N2 while r is about 30 per cent smaller and ^min 
is about 16 per cent larger representing a minimum mass ratio of 
about 1:4 rather than 1:5 for 03N2. There is a large overlap be¬ 
tween which galaxies are in the tail for N2 and for 03N2: 58 (76) 
per cent of the galaxies in the 03N2 (N2) tail are also present in the 
N2 ( 03N2) tail. _ 

iMaiolino et alJ ( 1200 Sl) find a relation of 12 -|- log(0/H) as 
a polynomial function of R23 = ([OIIJA3727 + [OIII]A4959 -F 
[OIII]A5007)/H/3. In this case we find a much wider and more 
skewed metallicity dispersion distribution with a scatter of a = 
0.80 dex and mean /r = 0.013 dex. The tail is bigger containing 
an excess of about 6.8 per cent of all galaxies. As with N2 there 
is a systematic offset in the metallicities of 0.05 dex but towards 
higher, rather than lower, metallicities. The high-SSFR peak of the 
galaxies in the tail is less prominent than for 03N2 causing the me¬ 
dian SSFR to be 0.08 dex lower (i.e. 0.16 dex lower than the median 
SSFR of the general population). The best-fitting values of the three 
model parameters are all larger than for 03N2, reflecting the wider 
tail, with a showing the biggest increase at 45 per cent. 70 (66) per 
cent of the galaxies in the 03N2 (R23) tail are also present in the 
R23 (03N2) tail. 

Finally, we look at the average of the N2 and R23 metallici¬ 
ties as was used in MIO. In this case, the scatter is slightly smaller 
at a = 0.045 dex while the mean is slightly larger at /r = 0.011 
dex. However, while the scatter is small the tail is large with an 
excess of 9.2 per cent of all galaxies. There is a systematic off¬ 
set towards higher metallicities of 0.26 dex. The SSFRs in the tail 
have the same bimodial shape as for 03N2 with an almost iden¬ 
tical median value as for 03N2 (0.01 dex lower than for 03N2, 
i.e. 0.09 dex lower than the general population). The best-fitting 
r and a are a bit higher than for 03N2 while ^min is significantly 
smaller representing a minimum mass ratio of about 1:8. Unlike 
the other calibrations for which the model can fit the residuals very 
well, in this case the fit is worsened by a few extreme outliers. 89 
(66) per cent of the galaxies in the 03N2 (N2-(-R23) tail are also 
present in the N2-I-R2 3 (03N2) tail. MIO used the N2 calib ration of 
IMaiolino et alj ( 1200 8h but we use the lPenicolo et alj (|2002|) calibra¬ 
tion because a sizeable fraction of the galaxies in o ur sample have 
1 ok([N II]A 6564/OIIA3727) < 1.2 for which the IMaiolino et al.l 
( 120081) calibratio n is not valid. However , for the subset of galaxies 
that do have valid iMaiolino et ^ ( l2008h N2 metallicities the resid¬ 


ual distribution is quite similar to the IMaiolino et al.l ( l2008h R23 
metallicities. 

From these comparisons, it is clear that while the shape of the 
metallicity dispersion around the FMR changes with the calibra¬ 
tion used, a low-r tail containing mostly the same low-metallicity 
galaxies is present across at least the three popular metallicity cali¬ 
brations based on 03N2, N2 and R23. Thus, our qualitative conclu¬ 
sion that the tail is consistent with reflecting the impact of mergers 
is unchanged and the uncertainties in the inferred model parameters 
arising from different metallicity calibrations seem to be no greater 
than the uncertainties caused by e.g. uncertainties in the merger rate 
or ignoring flybys. 


6.5 Aperture effects 


Spectroscopy is only available within the 3 arcsec aperture used 
for the SDSS fibre. Thus, we could only estimate metallicities 
for the central 3 arcsecs of the galaxies. The presence of metal¬ 
licity gradients could therefore potentially lead to many of the 
meta llicities being overes timated. To avoid significant aperture ef¬ 
fects iKewlev e t alj (l20 05h recommended a minimum redshift cut 
of a > 0.04 but Kewlev & Ellisonl ( l2008h found this to be insuffi¬ 
cient for galaxies with masses M* > However, for our 

sample we use a significantly higher redshift limit of 2 > 0.07. 
As the angular diameter distance at 2 = 0.07 is about 70 per cent 
greater than at 2 = 0.04, our lower redshift limit should be suffi¬ 
cient that aperture effects on the metallicity are relatively small for 
all masses. 

The total SFRs used in this paper are the sum of the 
SFRs within the apert ure, estimated using the method of 
ISrinchmann et al.l (l2004ll (these are based on several emission 
lines, but mainly on Ha), and the SFR outside the aperture, esti¬ 
mated from the photometry based on the assumption that the SFR 
has the same colour dependence outside of the aperture as within it. 
In contrast, MIO estimated the SFR only wit hin the aperture from 
the Hq line, as described in lKennicutd (Il998h . These SFRs are less 
uncertain than the total SFRs but only account for the star forma¬ 
tion within a central diameter of 4-11 kpc of the galaxies in the 
redshift range used in MIO and this paper. 

To examine the impact of using the aperture-corrected SFRs 
we fit FMRs using aperture-only SFRs, in addition to changing the 
metallicity calibrations, in table [T] These SFRs are derived in the 
same way as in MIO, correcting for dust using the Baimer decre¬ 
ment. The effect of aperture correcting the SFRs is typically to pro¬ 
duce a larger but also more concave tail leading to higher best¬ 
fitting values of r and a and lower values of ^min. As with different 
calibrations, the tail largely consists of the same galaxies when us¬ 
ing the same calibration and switching between aperture-corrected 
and aperture-only SFRs, with 75-90 per cent overlap. Overall, the 
effect of aperture correcting the SFRs is typically smaller than the 
effect of changing the metallicity calibration. 


6.6 Parametrization of the FMR 

In this paper, we have followed MIO and expressed the relation be¬ 
tween stellar mass, SFR and metallicity as a double quadratic. Usu¬ 
ally the name “fundamental metallicity relation” refers specifically 
to this form. While this is a widely used parametrization, the ques¬ 
tion of whether it is the most appropriate form is unresolved. Promi- 
nent alternatives includ e , for example, the “Fundamental Plane” 
dUara-Lopez et al.ll2010l . ItoT^) and the relation based on the an- 














































12 A. E. Gr0nnow, K. Finlator and L. Christensen 


alytical m odel of Lill y et al] l l2013h . The relation given by equa¬ 
tion 40 of iLillv et alj 1 2013h closely follows the MIO relation for 
the sample used in MIO by construction; it is only out s ide th e ob¬ 
served range of masses and SFRs that the yil^_et ^ (l2013h and 
MIO formulations differ significantly jMaier et~ .I l2014h . As our 
sample resemble s the MIO sample, differences with respect to the 

[Lilly et ^ ( 1201 3h relation will therefore be minimal. _ 

_ In the Fundamental Plane introduced in iLara-Lopez et al.l 

( l2010h the quadratic terms are dropped, leading to a linear rela - 
tion between mass, SFR and metallicity. iLara-Lopez et al.l l l2013l) 
argues that the flattening in the MIO relation at high masses is an 
artefact caused by the line ratios used in their metallicity calibration 
saturating at high metallicities. In this case, the underlying, physi¬ 
cal relation could well be a simple plane. We do not attempt to eval¬ 
uate whether the FMR is best describe d as a double quadratic or a 
plane (for a thorough comparison, see Ide los Reyes et al] (l2015h l. 
However, our merger model should also be able to fit the residu¬ 
als reasonably well when a plane is used. Fitting a plane to our 
stellar masses, SFRs and metallicities using least squares, we find 
that the distribution of residuals has the same shape as for the MIO 
relation, i.e., the residuals follow a Gaussian distribution with an 
overabundance of galaxies with lower metallicities than predicted. 
The tail appears to consist of largely the same population of galax¬ 
ies as when the MIO relation is used, i.e., the galaxies in the tail 
have slightly lower stellar masses, bimodial SFRs, lower metallici¬ 
ties and smaller half light radii than the general sample. 

Finally, we have also fitted an MZR of the form given in 
iTremonti et alj l l2004h . i.e. a fourth-order polynomial of the mass 
with no SFR dependence. In this case we also find that the residu¬ 
als follow the shape of a Gaussian distribution with a wing towards 
low metallicities. 

Generalizing the result that the tail could trace mergers to 
alternative functional dependences of metallicity on stellar mass, 
SFR and even further parameters is an avenue for future work. This 
approach will be particularly useful once complementary insight is 
available fro m theoretical models that treat both small fluctuations 
in gas flows dLillv et al.l20131 [Forbes et'^l2014l) and galaxy inter¬ 
actions. 


7 CONCLUSIONS 


Fitting FMRs to our large sample of SDSS galaxies, we have 
shown that a tail in the distribution of the FMR residuals towards 
lower metallicities is present regardless of the metallicity calibra¬ 
tion used. We develop a simple model where the tail consists of 
galaxies that are merging or have recently merged and show that 
this model is able to successfully reproduce the observed distri¬ 
bution of FMR residuals. The metallicity depressions and dilution 
time-scales that we derive by fitting this model are in good agree¬ 
ment with the results of merger simulations and observations of 
galaxy pairs. The success of this model suggests that the galaxies 
in the tail are recent mergers. This conclusion is supported by the 
fact that the galaxies in the tail form a distinct population with en¬ 
hanced star formation within the aperture where metallicity is also 
measured. 

We find that the average metallicity depression caused by a 
1:1 merger is abo ut 0.25 dex in agreeme nt w ith the hydrodynam - 
ical simulations of iMontuori et al. ( 1201 Ol) an d iRupke et alj j2010h. 
We also find the average depression of all 1:5 to 1:1 mergers to be 
0.114 dex which is consistent with the actual metallicity difference 
of the tail compared to the entire sample. This is a greater dilu¬ 


tion than what i s found in the metalli ci ty measurements of S DSS 
pair samples of Ellison et alJ (l2008bh . IScudder et al.l ( 1201 2h and 
iMichel-Eiansac et alj ( l200^r but our value is consistent with those 
studies when the differences between the minimum mass ratios in 
those samples and in our model are taken into account. 


We find that the average metallicity depression time-scale (the 
time from the onset of metallicity dilution until recovery to the pre¬ 
merger value) due to a m erger is 1.57 Gy r in go od agreement with 
the merger simulations of iMontuori et al.l (l2010l) . 


Currently the quantitative knowledge of the rate and effect of 
galaxy flybys is very poor. However in future as our understanding 
of flybys increases, the precision of our model can be significantly 
improved by extending the formalism to include flybys in addition 
to mergers. In addition, future hydrodynamical merger simulations 
more comprehensive than the merger simulations done so far, that 
probe the entire range of galaxy masses and mass ratios could allow 
for a more specific estimation of the dilution time-scale. Provided 
such a simulation is run for a sufficiently long time so that any 
potential stellar mass and mass ratio dependence of the metallicity 
depression time-scale can be found it would allow our model to pre¬ 
dict this time-scale for mergers of different member masses rather 
than as averaged over all mergers. In addition such a simulation 
might probe whether the mass ratio dependence of the metallicity 
depression (as time-averaged between pericentric passages) is ap¬ 
proximately linear as we assume or has a more complicated form. 


To further test if the galaxies in the tail are typically interact¬ 
ing, a morphological study of these galaxies could be undertaken, 
calculating merger indicators such as asymmetries or Gini coeffi¬ 
cients of the luminosity distributions. 


Due to the 3 arcsec aperture used in the SDSS, the accuracy 
of the galaxy parameters we use is limited by having to partly base 
stellar masses and SFRs on photometry, rather than spectroscopy 
only, and metallicities being estimated only within this aperture. 
In future, integral field spectroscopy should become available for a 
large fraction of our sample allowing us to obtain much more accu¬ 
rate global galaxy pa rameters. An impor tant step in this direction is 
the MaNGA survey dBundv et alj|2015]) , which is currently under¬ 
way and will obtain integral field spectroscopy of a total of 10,000 
nearby SDSS galaxies. While this only constitutes a relatively small 
fraction of our sample, it represents a factor ~ 100 increase in sam¬ 
ple size compared to previous integral field spectroscopy surveys. 


Our model makes some predictions that might be tested ob- 
servationally. Low-metallicity outliers from the FMR should have a 
higher pair fraction than the general population and the mass ratios 
of pairs should be closer to unity as one moves to greater metallic¬ 
ity depressions. On the theoretical side, large-scale hydrodynamical 
simulations that track the metallicity of thousands of g alaxies such 
as Illustris ( Genel et aTll2014l : IVogelsberger et al.l201^ or EAGLE 
dSchave et al.ll2015ll should, in principle, enable a statistical study 
of merger induced metallicity dilution that could test whether the 
galaxies in the low-metallicity tail are really interacting (or have 
recently interacted) and in that case compare their inferred dilu¬ 
tion time-scale, magnitude of metallicity depression and mass ratio 
dependence of those with our results. This is probably beyond the 
first generation of large-scale hydrodynamical simulations though, 
as outflows in Illustris are unphysic ally metal-depleted and EAGLE 
does not reproduce the MZR well dSchave et al.l201^ . 
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APPENDIX A: ESTIMATING THE POTENTIAL IMPACT 
OP SECONDARY MERGER MEMBERS 


Given the merger rate per galaxy per stellar mass ratio per lookback 
time for primary merger members, i.e. Si 1, 

and the distribution of galaxy stellar masses (j>{M,) (defined such 
that dM is the number of galaxies per volume with stel¬ 

lar mass between M and M + dM) we can calculate the rate of 
mergers per galaxy that galaxies of a given mass experience as the 
secondary member, i.e. > 1. 


mergers 
time • galaxy 


(^» > 1, Mo) = 


X 


Mq+SM 


()>(Mpri) f. 


Mq / 

(Mo+5M)/M 


d^P 

dt 


dMpr 


jMo+SM^^M)dM 


(Al) 


where the merger rate is evaluated at Mpu, and t = 1.30 Gyr 
which corresponds to « = 0.1, the typical reds hift of our sample. 
We use the mass function of iBaldrv et alj ilOld) . If we assume that 
the metallicity of these galaxies is affected in the same way as the 
primary galaxies we can get a rough estimate of the overall error 
9 that ignoring secondary galaxies introduces in the dilution time- 
scale T. We define 


7?pn(M) = 


mergers 


time • galaxy 


(Cmin < ^» < 1, M) 


and 


ftec(M) 


mergers 
time • galaxy 


(C* > hM) 


and calculate 9 by integrating the fraction of merging galaxies that 
are primary at each mass over the mass distribution of our sample. 


9 = 


'T’con-ected 


T 


J 0 d M Ppri + Pst 

f°° ^dM 
Jo dM 


■dM 


(A2) 


For tbe best-fitting value of ^min = 0.205 we find 9 — 0.29 so r 
would be overestimated by a factor ~ 3 — 4. 













